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We present new measurements of differential cross sections for Z/^*(-^ fifi) +jet+X production 
in a f fb _1 data sample collected with the DO detector in pp collisions at ^fs = f .96 TeV. Results 
include the first measurements differential in the Z/'y* transverse momentum and rapidity, as well 
as new measurements differential in the leading jet transverse momentum and rapidity. Next- 
to-leading order perturbative QCD predictions are compared to the measurements, and reasonable 
agreement is observed, except in the region of low Z/<y* transverse momentum. Predictions from two 
event generators based on matrix elements and parton showers, and one pure parton shower event 
generator are also compared to the measurements. These show significant overall normalization 
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differences to the data and have varied success in describing the shape of the distributions. 



PACS numbers: 12.38.Qk, 13.85.Qk, 13.87.-a 

The production of W or Z bosons in association with 
jets is an important signal at hadron colliders such as 
the Fermilab Tevatron Collider and the CERN Large 
Hadron Collider. Leptonic boson decays can be iden- 
tified with little background, and measurements of the 
boson and jet kinematics provide good tests of pertur- 
bative QCD (pQCD) calculations and modeling. Such 
events also form the main background to many processes 
with much smaller cross sections, including production 
of the top quark, Higgs boson, and particles expected 
in some supersymmetric scenarios. Accurate theoretical 
modeling of W or Z boson + jet final states is a key 
element in studying such rare processes; developing and 
testing these models relies upon input from experimental 
measurements of boson + jet production. 

Previous boson + jet measurements at the Fermi- 
lab Tevatron Collider [D, 0] have included comparisons 
with next-to-leading order (NLO) pQCD predictions 
from MCFM 0] and a tree-level matrix element calcula- 
tion matched to a parton-shower (ME+PS) Monte Carlo 
event generator [4J . In this Letter we describe new mea- 
surements of differential cross sections in Z/j*(— > [i\x) + 
jet + X production in pp collisions at y/s = 1.96 TeV, 
with a data sample corresponding to 0.97 ± 0.06 fb _1 [5] 
recorded by the DO detector between April 2002 and 
February 2006. The measurement is carried out in 
a region of dimuon mass 65 < M MAt < 115 GeV in 
which the inclusive cross section for Z/j* production 
is approximately equal to that of pure Z boson pro- 
duction, and the measured distributions are corrected 
to the particle level [6J]. Differential cross sections 
binned in the Z/j* momentum component perpendicu- 
lar (transverse) to the beam, p%, (daz+jet+x I&Pt) an d 
rapidity Q, y z (daz+jct+x /dy z ), are presented here 
for the first time, and binned in the leading (in px) 
jet pt (daz+jet+x/d-Pr') extended to lower p^ , and y 
(dcrz+jot+x/dy jct ) covering a wider y jct range than pre- 
vious studies. Comparisons are made using NLO pQCD 
predictions from MCFM with non-perturbative corrections 
applied. Currently the best tools for generating simu- 
lated boson + jets events are tree-level ME+PS calcula- 
tions. Several such calculations are available using differ- 
ent schemes to combine the matrix element and partons 
shower contributions, and with some differences in the 
predicted kinematics Q. Two such ME+PS event gener- 
ators are com par ed to the measurements: ALPGEN Q, 
with PYTHIA [lfj used for the parton-showering; and 
SHERPA [llj]. Finally, a pure parton-shower prediction 
from pythia is also compared. 

The measurements are made with the DO detector, 
which is described in detail elsewhere a brief 



overview is given here of the most relevant components 
for this analysis. The interaction region is surrounded by 
a magnetic tracking system, comprising a silicon micro- 
strip tracker and a fiber tracker, both located within 
a 2 T superconducting solenoid magnet. Three liquid- 
argon/uranium calorimeters surround the tracking sys- 
tem: a central section covering pseudo-rapidity \r]\ < 
1.1 [l3j|, and two end calorimeters that extend coverage 
to 1.1 < |r/| < 4.2, each housed in separate cryostats. 
Scintillators between the cryostats sample shower energy 
for 1.1 < |r/| < 1.4. Luminosity is measured using plastic 
scintillator arrays located in front of the end calorime- 
ter cryostats, covering 2.7 < |?7| < 4.4. A muon sys- 
tem surrounds the calorimetry, consisting of three lay- 
ers of tracking detectors and scintillation trigger coun- 
ters; these provide muon identification and triggering for 
1 77 1 < 2. Sandwiched between the first and second layer 
are 1.8 T toroidal iron magnets, allowing an independent 
momentum measurement in the muon system. 

Events used in this measurement are selected with a 
suite of triggers using information from the muon and 
tracking systems and are required to have two muon can- 
didates reconstructed in those systems. The primary 
collision vertex in each event is reconstructed requir- 
ing at least three tracks and applying fit quality cuts. 
To reject mis-reconstructed events and cosmic rays, the 
muons must be consistent with this vertex in directions 
both transverse and parallel to the beam. Based on 
the information from the tracking system, the muons 
are required to have pr > 15 GeV and dimuon mass 
65 < < 115 GeV. 

Jets are reconstructed from clusters of energy de- 
posited in calorimeter cells using the DO Run II midpoint 
cone algorithm [14| with a splitting/merging fraction of 
0.5 and a cone size of ATI = y/(A<j)) 2 + (Ay) 2 = 0.5, 
where </> is the azimuthal angle; jets caused by noise 
are rejected with quality and shape cuts. Jets are cor- 
rected for the calorimeter response, instrumental out-of- 
cone showering effects, and additional energy deposits 
in the calorimeter that arise through detector noise and 
pile-up from multiple interactions and previous pp bunch 
crossings. These jet energy scale corrections are deter- 
mined using transverse momentum imbalance in 7 + jet 
events, after the electromagnetic response is calibrated 
using Z/7*— > ee events. For clarity in the following dis- 
cussions, the measured jet transverse momentum and ra- 
pidity after these corrections are denoted PjP T and y JET , 
to distinguish them from the particle level quantities Pj, 
and j/ jct . 

Further selections ensure that the measurement is car- 
ried out in regions with high acceptance and well un- 



5 



derstood detector performance: the muons are required 
to have \rj\ < 1.7, the primary vertex must be within 
50 cm of the detector center along the beam direction, 
and jets are required to have a p^ ET > 20 GeV and 
|yJET| < 2_g, Additionally, events with jets in the p T ET 
range 15-20 GeV are kept in the sample for studies of 
the effects of detector resolution. 

The main source of background in this analysis is 
muons from semi-leptonic decays in high energy jets or 
W+jet production. This is reduced to negligible lev- 
els (< 0.5% of the final sample) by limiting the sum 
of track momenta and the calorimeter energy allowed 
around each muon. The muons are also required to 
not overlap with any jet by requiring angular separation 
v / (A0(^i,jet)) 2 + (Ar](fi, jet)) 2 > 0.5. Other sources of 
background (e.g., top quark production, Z/j* — ► t + t~) 
are estimated using simulation and found to be negligi- 
ble (< 0.1%). A total of 59336 Z/j* /x+^t" candidate 
events are selected before jet requirements, of which 9927 
contain at least one jet with PjP T > 20 GeV passing all 
selections. 

Two simulations of Z/-y*+jets events are used: a sam- 
ple generated with PYTHIA v6.323, and a sample gener- 
ated with ALPGEN v2.05 using PYTHIA for parton shower- 
ing, both with the PYTHIA underlying event model config- 
ured usin g tu ne A [HI . Both samples are passed through 
a GEANT [16( simulation of the detector response. Real 
data events from random bunch crossings are overlaid 
on the simulation to reproduce the effects of multiple pp 
interactions and detector noise. 

The muon trigger efficiency is measured in data and pa- 
rameterized in terms of variables related to the geometry 
of the muon and tracking systems. The selected events 
are then corrected on an event-by-event basis for this ef- 
ficiency, with the average efficiency being (88.3 ± 0.3)%, 
quoting just the statistical uncertainty. Muon recon- 
struction, tracking, and isolation efficiencies are mea- 
sured in the data and in the simulation; scale factors are 
applied to correct for the differences. The total system- 
atic uncertainty on the muon trigger and identification 
efficiency translates into a 5% uncertainty on the mea- 
sured cross sections, with no significant dependence on 
the variables studied in this analysis. Transverse mo- 
mentum imbalance in Z/j*— > ee + jet events is studied 
in data and simulation, and factors applied to the simu- 
lation to correct the jet response for any differences. 

To extract the differential cross sections, we correct 
the reconstructed data distributions to particle-level dis- 
tributions, deriving the corrections from simulation. We 
first select Z/j* plus jet events based on the detector 
response in simulated events. These Z/j* and jet vari- 
ables are compared to quantities independently measured 
directly from the particles in the simulated events, apply- 
ing comparable kinematic selections to minimize accep- 
tance effects. For this particle level selection, the Z/j* 
mass and kinematics are reconstructed from the gener- 



ated muons after QED final state radiation (FSR), re- 
quiring 1 77^ | < 1.7 and 65 < Af w < 115 GeV. Jets are 
reconstructed using the same reconstruction algorithm, 
now on all final state particles excluding the Zj 7* decay 
products, and selected requiring |y jet | < 2.8. These par- 
ticle jets are matched to jets reconstructed in the sim- 
ulation by requiring A7£ < 0.5. We quote results for 
leading particle jets with p T :t > 20 GeV; however, due to 
instrumental effects and resolution, measured jets with 
p T ET > 20 GeV include significant contributions from par- 
ticle jets with lower p!^. To study this effect, jets at the 
particle level are reconstructed to very low px (3 GeV) . 

We next describe the process of correcting the px dis- 
tribution of the leading (in px) jet from the measured 
level to the particle level; the treatment of pj, is very 
similar. The main complexity in the jet px corrections 
arises from the experimental px resolution affecting the 
relationship between particle jets and the corresponding 
jets reconstructed in the detector. First, the finite energy 
resolution can change the px ordering of jets between the 
particle level and detector level. To account for this we 
correct the measured p T ET distribution to remove lead- 
ing measured jets matched to sub-leading particle jets, 
based on a study of simulated events. Here we also re- 
move measured jets arising from additional collisions in 
the event, modeled by the random bunch crossings from 
real data overlaid on the simulation. This combined cor- 
rection averages (11.8±0.2)%, varying from (33.8±0.2)% 
in the range 15 < p T ET < 20 GeV, to (3.1 ± 0.1)% for 
p T ET > 50 GeV. The second effect of the resolution re- 
sults in some jets from a given p^ bin being measured 
in a different p T ET bin. This effect is mitigated to a de- 
gree by the choice of binning for the measurement: bins 
are taken to be wider than the detector resolution and 
to contain a sufficient number of events so that statisti- 
cal fluctuations do not dominate the final uncertainty on 
each bin. Studying the p^? 1 and the corresponding mea- 
sured p T ET f° r jets in the full detector simulation allows 
the remaining effect to be parameterized in a "migration 
matrix" (see Fig. [T]), with element i,j being the prob- 
ability for a particle jet in p!^ bin i to be measured in 
p T ET bin j. The data distribution is then corrected using 
a regularized inversion of this matrix 17J , with the con- 
straint that the resulting distribution does not have large 
second derivatives. Including the reconstructed jets with 
15 < p T ET < 20 GeV in the matrix further constrains 
the effects of lower p 3 ^ particle jets fluctuating up in re- 
constructed p T ET ■ Finally, the distribution is corrected 
for efficiency and acceptance calculated from simulation, 
then divided by the bin widths and integrated luminosity 
to yield the differential cross section. 

Uncertainties on the differential cross section are de- 
rived empirically through ensemble tests. A set of 100 
ensembles of the same size as the data set are drawn 
from a pythia sample. To reproduce the p!^. distribu- 
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FIG. 1: The migration matrix for leading pi^ . Element i,j is 
the probability for a particle jet in pr bin i to be measured 
in pr bin j, represented by the area of each box. Each row 
sums to unity. 



tion in data, the PYTHIA p^ spectrum is re-weighted us- 
ing a function derived from the corrected data and a large 
(2.5 million events) independent PYTHIA sample. Apply- 
ing this function to the ensembles reproduces the data 
Pj, distribution while retaining realistic statistical fluc- 
tuations. The measured distribution in each ensemble 
is then corrected in the same way as the data. Uncer- 
tainties are extracted by taking the fractional difference 
between the fully corrected distribution and the actual 
Pj? distribution in each ensemble. The systematic un- 
certainty is the mean fractional difference in each p!£ 
bin over all 100 ensembles; the statistical uncertainty is 
the RMS around the mean. These statistical uncertain- 
ties account for the statistics in each measured bin, and 
the effects on those statistics of migrations between bins. 
The systematic uncertainties are typically below 2%, and 
the statistical uncertainty in each bin varies from 2% at 
low to 11% at high p^*. 

Further systematic uncertainties are then assessed. 
Varying the re- weighting function used in generating the 
ensembles produces uncertainties at the 3% level, mostly 
at low p^* due to the weaker constraints on the particle 
jet spectrum below the measured region. Studies of the 
jet resolution and reconstruction efficiency show small 
effects, and larger (< 3%) effects are seen by varying 
the jet energy scale within uncertainties in the data and 
simulation. All other systematic uncertainties studied 
produced negligible effects. No strong correlations are 
observed between the various sources of uncertainty, and 
the individual contributions are combined in quadrature 



to obtain the total systematic uncertainty. 

The p^ distribution is corrected using the same ap- 
proach, employing a regularized inversion of the migra- 
tion matrix with uncertainties derived from ensemble 
testing; in this case, the muon pt resolution is the source 
of migration. Along with the sources of systematic uncer- 
tainty considered for the leading jet pr, the uncertainty 
on the agreement between the muon resolution in simu- 
lation and data is also considered. Varying the resolution 
in simulation within uncertainties produces effects below 
(2-3)% on the differential cross section. Varying the jet 
energy scale produces systematic effects of up to 10% in 
the region of pf < 20 GeV, which is sensitive to jets close 
to the reconstructed p T ET cutoff. 

The measurements of y z and y^ ct are significantly less 
challenging, and the method used on these variables is 
covered briefly. These distributions do not suffer from 
significant resolution effects on the rapidity measure- 
ment, but still need to be corrected for efficiency and 
acceptance. To do this, the ratio of particle level to mea- 
sured events in each bin is calculated in simulation and 
applied to the measured data distribution. Ensemble 
testing is then used to measure the uncertainties, with 
the same sources as the p J ^* and p^ measurements re- 
spectively. Including the jet systematics covers the cor- 
relations between the rapidity and Pj. distributions, tak- 
ing into account changes in the rapidity distributions as 
events enter or leave the sample due to jet migrations 
across the p^ ET selection of 20 GeV. As the distributions 
are symmetric around zero, \y\ is measured in both cases 
to increase the statistics in each bin. 

Integrating over any of the differential cross sections 
yields the ZJ^*{— * /x/i)+jet+X cross section, which we 
measure to be 18.7 ± 0.2(stat.) ± 0.8(syst.) ± 0.9(muon) ± 
l.l(lumi.) pb, with the following requirements: all bo- 
son properties are calculated from the muons after QED 
FSR, and the muons are required to have \y\ < 1.7 and 
dimuon mass 65 < M w < 115 GeV; particle jets are 
reconstructed using the DO Run II midpoint algorithm 
with a splitting/merging fraction of 0.5 and a cone size 
of A1Z = 0.5 on all final state particles except the Z/j* 
decay products and any FSR photons from the muons, 
and are required to have |y jet | < 2.8 and Pj> t > 20 GeV. 
The quoted muon uncertainty covers the muon identi- 
fication and trigger efficiency determination. Different 
definitions of observables complicate comparisons, but 
this represents a significant reduction in cross section 
uncertainty from the previously published DO result us- 
ing 0.4 fb _1 of data [l(, and is of comparable accuracy 
to the CDF result using 1.7 fb _1 2]. For reference, 
we also measure an inclusive Z/j*(— > p,p) cross sec- 
tion requiring only 65 < M w < 115 GeV, taking the 
muons after QED FSR with no rapidity requirements, 
and no jet requirements, to be 233 ± l(stat.) ± 8(syst.) ± 
12(muon) ± 14(lumi.) pb. Adding only the requirement 
that the muons have |y| < 1.7 yields a cross section of 
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TABLE I: Non-perturbative and FSR correction factors ap- 
plied to the mcfm prediction for leading p^* in Z/7* + jet + 
X events. 



jet 

P T 


Non-pert. 


FSR 


(GeV) 


Corr. 


Corr. 


20-30 


1.041 


0.977 


30-40 


1.017 


0.977 


40-50 


1.001 


0.977 


50-60 


0.995 


0.977 


60-70 


0.991 


0.977 


70-80 


0.989 


0.978 


80 - 100 


0.988 


0.978 


100-130 


0.986 


0.978 


130 - 200 


0.984 


0.978 



118 ± 0.5(stat.) ± 4(syst.) ± 6(muon) ± 7(lumi.) pb. The 
muon and luminosity uncertainties are completely corre- 
lated between the inclusive Z/j* (with and without the 
muon y requirement) and Z/7*+jet measurements, and 
all other systematics are uncorrelated. 

A number of theoretical predictions are now compared 
to the measured integrated and differential Z/-f*(-^ 
Hii)+\et+X cross sections. NLO pQCD calculations are 
obtained using mcfm together with the NLO CTEQ6.6M 
parton distribution functions of the proton (PDF) [lil ]. 
The associated PDF error sets are used to assess un- 
certainties, which are about 3%. Renormalization and 
factorization scales are set to the sum in quadrature of 
the mass and pt of the Z boson, and uncertainties are 
derived by varying both scales down or up together by 
a factor of two, which changes the prediction by ±7%. 
Non-perturbative corrections for hadronization and the 
underlying event are derived from pythia v6.418 Z/j* 
+ jet production, with the leading order (LO) PDF 
CTEQ5L and the underlying event tune DW [H. 
These are derived by comparing the full prediction (taken 
from the final state particles, including the underlying 
event) to the purely perturbative part (calculated from 
partons taken after the parton shower, with no underly- 
ing event). Corrections for QED FSR from the muons are 
derived from the same PYTHIA sample, by comparing the 
prediction calculated using the muons after FSR to those 
using the muons before FSR. These FSR corrections are 
around 2% caused mainly by events migrating out of 
the mass window, with little dependence on any vari- 
able considered except at low and high y z . The non- 
perturbative and FSR corrections are given in Tables [H 
HTlHTll and[lV] The prediction for the Z/V+jet+A" cross 
section from NLO pQCD with our stated acceptance cuts 
and after corrections (hereafter referred to as the NLO 
pQCD prediction), is 17.3 ± 1.2(scale) ± 0.5(PDF) pb. 
This is 5% below the measured value, and within uncer- 
tainties. For reference, mcfm is also used to calculate the 
LO pQCD prediction. After acceptance cuts and correc- 
tions this is 12. (scale) ± 0.3(PDF) pb. 



TABLE II: Non-perturbative and FSR correction factors ap- 
plied to the mcfm prediction for leading \y sat \ in Z/y* + jet 
+ X events. 



\y 


l\Tr\n t~i ot'L 
1>I UI1 JJcI L . 


-T on 




Corr. 


Corr. 


0.0-0.4 


1.026 


0.977 


0.4-0.8 


1.025 


0.977 


0.8-1.2 


1.022 


0.977 


1.2-1.6 


1.018 


0.977 


1.6-2.0 


1.001 


0.977 


2.0-2.4 


1.009 


0.977 


2.4-2.8 


0.983 


0.977 


TABLE III: Non-perturbative and FSR correction factors ap- 


plied to the mcfm prediction for p z in Z/7* + jet + X events. 


Pt 


Non-pert. 


FSR 


(GeV) 


Corr. 


Corr. 


0-10 


1.099 


1.268 


10-18 


1.284 


1.041 


18-26 


1.021 


0.977 


26-35 


0.995 


0.972 


35-45 


0.997 


0.973 


45-60 


0.999 


0.967 


60-80 


1.000 


0.959 


80 - 120 


1.000 


0.951 


120-200 


1.000 


0.944 


Comparisons are also made using 


three event genera- 


tors: i) a sample 


generated with ALPGEN v2.13 with up 


to three partons 


in the matrix element calculation, and 


the factorization and renormalization scales squared set 


to the sum in quadrature of the mass and pr of the Z bo- 


son, and pythia with tune QW [15| 


used for the parton 


showering; ii) a 


sample generated with SHERPA vl.1.1, 


again with up to three partons in 


the matrix clement 


calculation, and the default parton showering algorithm 


(apacic) and underlying event model (amisic); iii) an in- 


clusive Z/j* sample generated with pythia v6.418, using 


tune QW. In order to use consistent PDF with all three, 


TABLE IV: Non-perturbative and FSR correction factors ap- 


plied to the mcfm prediction for \y z \ 


in Z/7* + jet + X 


events. 






\y z \ 


Non-pert. 


FSR 




Corr. 


Corr. 


0.0-0.2 


1.024 


0.979 


0.2-0.4 


1.024 


0.978 


0.4-0.6 


1.024 


0.978 


0.6-0.8 


1.021 


0.978 


0.8-1.0 


1.017 


0.977 


1.0-1.2 


1.013 


0.975 


1.2-1.4 


1.010 


0.971 


1.4-1.6 


1.001 


0.965 


1.6-1.8 


1.000 


0.928 
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the NLO CTEQ6.1M [20| is chosen. For both ALPGEN 
and SHERPA, jets from the matrix element calculation 
are required to have px > 15 GcV, and a separation 
AB, > 0.4, to ensure full coverage of the measured phase 
space. For all three event generators, the boson kine- 
matics are calculated from the muons after QED FSR, 
for consistency with the measured observables. After ap- 
plying our stated particle level Z/-f*+ jet selection, the 
predicted cross sections are 11.6 pb (alpgen), 15.0 pb 
(sherpa), and 12.1 pb (pythia). 



D0 Run II, L=1.0fb 1 
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FIG. 2: (a) The measured cross section in bins of leading 
for Z/'j* + jet + X events. Predictions from NLO pQCD 
and ALPGEN are compared to the data, (b) The ratio of data 
and predictions from NLO pQCD + corrections, sherpa, and 
pythia to the prediction from ALPGEN. 

The differential cross sections are shown binned in 
leading pij. (Fig. [5]), leading jet y (Fig. [3]), in events 
with at least one jet (Fig. and y z in events with 
at least one jet (Fig. [5]). Data points in each bin are 
placed where the differential cross section in simulation 
is equal to the bin average [21]. The data are shown 
with statistical uncertainties (inner error bar) and sum 
in quadrature of statistical and systematic uncertainties 
(outer error bar) , excluding the uncertainties on the mea- 
sured integrated luminosity and the muon identification 
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FIG. 3: (a) The measured cross section in bins of leading \y' ct \ 
for Z/'y* + jet + X events. Predictions from NLO pQCD 
and ALPGEN are compared to the data, (b) The ratio of data 
and predictions from NLO pQCD + corrections, sherpa, and 
pythia to the prediction from ALPGEN. 



and trigger efficiencies. These final two uncertainties are 
completely correlated between bins and with the muon 
and luminosity uncertainties on the measured inclusive 
Z l^j* cross section; however, they are included to form 
the total uncertainty, shown as the shaded region. For 
clarity, only the predictions of NLO pQCD and ALPGEN 
are shown in part (a) of each figure, though the prediction 
from NLO pQCD is not shown at low (Fig. 2]) where 
non-perturbative processes dominate over the NLO con- 
tribution. The data results are also provided in Tables fVl 
VH EED and [ylTTI In part (b) of each figure, the distri- 
butions from data, NLO pQCD, sherpa and pythia are 
shown divided by the prediction from alpgen. The NLO 
pQCD prediction is shown with the scale and PDF uncer- 
tainties combined in quadrature as a hatched region; the 
scale uncertainty is approximately a factor of two larger 
than the PDF uncertainty across all distributions. 

In summary, we have measured differential cross sec- 
tions for Z/7*+jet+X production with 0.97 ± 0.06 ftr 1 
of integrated luminosity recorded by the DO experiment 
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FIG. 4: (a) The measured cross section in bins of for 
Z/'y* + jet + X events. Predictions from NLO pQCD and 
ALPGEN are compared to the data, (b) The ratio of data 
and predictions from NLO pQCD + corrections, SHERPA, and 
pythia to the prediction from ALPGEN. 



FIG. 5: (a) The measured cross section in bins of \y z \ for 
Z/'y* + jet + X events. Predictions from NLO pQCD and 
ALPGEN are compared to the data, (b) The ratio of data 
and predictions from NLO pQCD + corrections, SHERPA, and 
pythia to the prediction from ALPGEN. 



in pp collisions at y/s = 1.96 TeV. We presented the first 
results binned inp^ and y z ; as well as new results binned 
in leading jet pr and y, extending the measured pj? and 
yi ct ranges substantially. The total Z/^+jet+X cross 
section measured in data is 5% above the prediction from 
NLO pQCD + corrections, which is within the total un- 
certainties. This is comparable with the trend observed 
in previous measurements Q], although direct compar- 
isons are complicated by different definitions of the final 
observables. The shapes of the differential distributions 
are generally well described by NLO pQCD, though the 
distribution at lower pj, (below the cutoff of 20 GeV) 
is dominated by non-perturbative processes. The total 
cross sections predicted by the ALPGEN and pythia event 
generators are significantly below the measured value and 
are more consistent with the LO pQCD predictions. The 
prediction from SHERPA lies between the LO and NLO 
pQCD prediction. The shapes of the data distributions 
are generally well described by ALPGEN, except at low 
p^. There is also indication that the jet rapidity distri- 



bution is narrower in ALPGEN than in data, NLO pQCD, 
SHERPA and pythia. Comparisons to the other event 
generators show that SHERPA has a slope in and p z 
relative to the data, with more events at high px com- 
pared to low pt; pythia shows the opposite behavior. 
This measurement tests the current best predictions for 
heavy boson + jet production at hadron colliders. As 
the data are fully corrected for instrumental effects, they 
can be directly used in testing and improving the existing 
event generators, or any future calculations and models. 
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TABLE V: The measured cross section in bins of leading 
for Z/-y* + jet + X events. Uncertainties are split into statis- 
tical and systematic; these are combined with an additional 
constant 8.0% normalization uncertainty from the luminosity, 
trigger, and muon identification to form the total uncertainty. 



Pt 


Bin ctr. 


da /dp^ t 


oa sta t. 


OfT S y S t. 


0<T total 


(GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 


(%) 


20-30 


24. 7 


0.867 


2.4 


4.7 


9.6 


30-40 


34.8 


0.402 


2.7 


2.9 


8.9 


40-50 


44.8 


0.219 


3.2 


4.2 


9.5 


50-60 


54.8 


0.134 


4.4 


5.1 


10.5 


60-70 


64.8 


. 0854 


5.1 


4.6 


10.6 


70-80 


74.7 


. 0535 


6.3 


6.7 


12.2 


80 - 100 


89.0 


. 0301 


7.9 


6.2 


12.9 


100 - 130 


113.5 


0.0118 


8.4 


10.6 


15.7 


130 - 200 


157.7 


. 00226 


13.3 


18. 1 


23.8 



TABLE VI: The measured cross section in bins of leading | j/ Jct | 
for Z/-y* + jet + X events. Uncertainties are split into statis- 
tical and systematic; these are combined with an additional 
constant 8.0% normalization uncertainty from the luminosity, 
trigger, and muon identification to form the total uncertainty. 



1 „ jet 1 

\y \ 


Bin 


da/d|j/ jet | 


<5<7stat. 


<5<7 S yst. 


So-total 




center 


(Pb) 


(%) 


(%) 


(%) 


0.0-0.4 


0. 189 


12.01 


2.3 


2.9 


8.8 


0.4-0.8 


0.609 


10.66 


2.2 


2.9 


8.8 


0.8-1.2 


1.00 


9.13 


2.8 


5.3 


10.0 


1.2-1.6 


1 .40 


6.70 


3.2 


6.9 


11.0 


1.6-2.0 


1 .80 


4.38 


4.2 


7.6 


11.8 


2.0-2.4 


2.20 


2.46 


4.9 


11.1 


14.5 


2.4-2.8 


2.59 


1.40 


7.3 


13.1 


17.0 
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TABLE VII: The measured cross section in bins of for Z/-y* 
+ jet + X events. Uncertainties are split into statistical and 
systematic; these are combined with an additional constant 
8.0% normalization uncertainty from the luminosity, trigger, 
and muon identification to form the total uncertainty. 



Pt 


Bin ctr. 


do/dp? 


(Jfatat. 


<5o"syst. 


^CTtotal 


(GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 


(%) 


0-10 


5.2 


. 0410 


5.6 


18.5 


20. 1 


10-18 


14.5 


0. 151 


4.4 


15.0 


17.0 


18-26 


21.7 


0.448 


2.5 


9.5 


12.6 


26-35 


31.5 


0.525 


2.5 


6.8 


10.8 


35-45 


39.8 


0.342 


2.3 


2.2 


8.6 


45-60 


52.1 


0.179 


2.9 


4.9 


9.8 


60-80 


69.3 


. 0748 


3.7 


4.6 


9.9 


80-120 


97.3 


. 0233 


5.8 


3.0 


10.3 


120-200 


148.6 


. 00309 


10.8 


6.7 


15.0 



TABLE VIII: The measured cross section in bins of \y z \ for 
Z/'y* + jet + X events. Uncertainties are split into statis- 
tical and systematic; these are combined with an additional 
constant 8.0% normalization uncertainty from the luminosity, 
trigger, and muon identification to form the total uncertainty. 



\y z \ 


Bin 


do/d\y z \ 




5cr syst . 


^CTtotal 
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(Pb) 


(%) 


(%) 


(%) 
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6.5 


10.6 
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17.33 
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7.1 
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